Probing the susceptibility x'(q) in cuprates using Ni impurities 
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Recent NMR experiments on YBa 2 (Cui_ a; NL J )306 + ,5 (Bobroft et al. Phys. Rev. Lett. 79, 2117 
(1997)) have shown that Ni impurities provide an important probe for the momentum dependence 
of the static spin susceptibility Combining the measurements of the 17 line width Au with 

those of the Gaussian relaxation time T 2 g , we find that the magnetic correlation length £ must have 
a substantial temperature dependence. Using numerical simulations we present a detailed analysis 
of Av as a function of temperature and impurity concentration. For a Lorentzian shape of x'(q) Az/ 
strongly depends on £, in contrast to the Gaussian form. This result together with the experimental 
finding that £ is temperature dependent rules out the Gaussian form of x'(<3)- 
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Precise knowledge of the spin susceptibility x(<L w ) of 
the cuprates is essential for understanding their unusual 
normal state properties. The imaginary part, }("{(\,uj) 
can be probed either by inelastic neutron scattering (INS) 
[p-p], or in the low frequency limit by NMR measure- 
ments of the spin-lattice relaxation rate 1/Ti J5j. In con- 
trast, one knows little about the real part of the suscep- 
tibility, x'Cl)? since information can, so far, only be ex- 
tracted from the NMR observation of the Gaussian com- 
ponent of the transverse relaxation time, T2G, of planar 
Cu ||]7|]. In particular, the analysis of INS and NMR ex- 
periments has not yet led to a consensus on the shape of 
x(q, w) in momentum space and the temperature (T) de- 
pendence of the antiferromagnetic correlation length, £. 
In this communication we present new insight into this 
issue based on experiments by Bobroff et al. Our 
principal conclusions are that £ in YBa2Cu30g+5 is In- 
dependent and that the Lorentzian form of x' (<j) provides 
a completely consistent description of the data, whereas 
the Gaussian form can be ruled out. 

Bobroff et al. ||] recently presented a novel approach 
to the measurement of x'(q) using Ni impurities in 
YBa2(Cui_ ;c Ni x )306+<5. These impurities induce a spin 
polarization at the planar Cu sites via x'(q). The hy- 
pcrfinc coupling between Cu and O induces a spatially 
varying polarization and an additional broadening 



Az*, 



Ai/ - Av = af(0/T 



(1) 



of the planar 17 O NMR, where Av and A^o are the to- 
tal and x = line width, respectively. In Eq. [j] a is 
the overall amplitude of %'(q) and /(£) characterizes the 
dependency of Av on £ (a = 47rx* in the notation of 
Ref. (§ and 0). Finally, the factor 1/T is caused by 
the Curie behavior of the Ni impurities in YBa2Cu306+s 
P,fi~i| with effective moment p e s ~ 1.9/xb (1.59/is) for 
6 = 0.6(<5 = 1). Bobroff et al. found that TAv(T) 
strongly depends on temperature and the Ni concentra- 
tion x in the sample. Furthermore they observed a much 
stronger broadening in the undcrdoped, 8 = 0.6, sample 
than in the overdoped one with 5 = 1. Performing nu- 
merical simulations of the NMR line shape by assuming a 



Gaussian form for x'Cq)) they found that /(£) is basically 
constant for all physically reasonable values of £. Com- 
bining these results with T 2 g data by Takigawa JFJj , they 
concluded that £ is T-independent for the underdoped 
samples. On the other hand, in every scenario of cuprate 
superconductors in which the anomalous low-energy be- 
havior is driven by spin fluctuations one would expect 
the correlation length £ to be T-dependent (for recent 
reviews, see: JlJ,|l^|). Thus their result has important 
implications about the mechanism of superconductivity. 
We recently pointed out |[o]], that our simulations using 
a Lorentzian form of x'(*l) yield a different result and are 
actually compatible with a T-dependent £. 

Before going into the details of our calculations, it 
is important to notice that the fact that £ must be In- 
dependent can be deduced even without a detailed model 
from the very experimental data by Bobroff et al. || for 
Az/(T) and Takigawa H] for T 2 g- To show this, we need 
to recognize that we can always express T 2 g as a product 
of a and a function of £, namely 



T, 
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We can then eliminate a by forming the product 

/(0 



TAv imp T2G — 
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(2) 
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which depends solely on £. In Fig. Jn we plot the product 
T 2 G? 1 A^i m p as a function of T p5| . We see that this 
product is strongly T-dependent, dropping by more than 
a factor of 2 between 100 K and 200 K . Therefore £ must 
have a substantial T-dependence. 

To have a more quantitative insight into the T- 
dependence of TAv(T) of Ref. ||, we must go into de- 
tails. We present in the following a theoretical analysis 
of the 1 O line shape using a method first applied by 
Bobroff et al., to simulate their experimental data. 

To simulate the 17 line shape numerically, we dis- 
tribute Ni impurities on a (100 x 100) lattice with concen- 
tration ix randomly at positions rj on a two dimensional 



lattice [M. We consider the Ni impurities as foreign 
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FIG. 1. TAi/i mp T2G as a function of temperature, us- 
ing the original (empty circles) or Ti corrected (filled cir- 
cles) T2G data |15| ]. This product is independent of a and 
its T-dependence must be caused by a T-dependent £. We 
used the data from Ref. [H and Q . 



atoms embedded in the pure material, which is charac- 
terized by a non-local spin-susceptibility x'(q)- I n the 
following Si characterizes the spin dynamics of the pure 
material, Si the difference at the Ni site brought about 
by the Ni. These Ni spins polarize the spin Sj of the 
itinerant strongly correlated electrons. To calculate the 
induced moments we need to know how the Ni impurities 
couple to these spins. Without discussing the microscopic 
origin of the effective Ni spin Sj , we assume that it obeys 
a Curie law, and that the coupling to the spin Sj occurs 
via an on-site interaction described by 



Hi 
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(4) 



The coupling constant J is an unknown parameter of the 
theory and will be estimated below. Furthermore, we will 
assume like Bobroff et al. that the Ni impurities do not 
change the magnetic correlation length or the magnitude 
of the spin susceptibility. 

For the NMR experiments we consider an exter- 
nal magnetic field Bq along the z-direction. The 
Ni spins have a non-zero average value obeying 



{SI) 



Ccurie-Bo/T with Curie constant Cc U rie 



Peff/(2v / 3fc B ) 10. 

Adopting a mean field picture, the induced polariza- 
tion for the electron spins at the Cu sites is given by 



J 
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Here, x'(r) is the real space Fourier transform of x'(q). In 
the following we consider two different forms of the spin 
susceptibility |0. For the commensurate case, there is 



only one peak, whereas in the incommensurate case, one 
has to sum over four peaks. The Gaussian form of x'(q) 
is given by 



X^(q)=< 2 exp(-(q-Q) 2 £ 2 ) 
and the Lorentzian form by 

xL(q) = ^ 2 /(i + (q-Q) 2 a- 



(6) 
(7) 



Since the question whether there exist incommensurate 
peaks in YBa2Cu30g+5 has not been settled yet 0], we 
will consider below both cases, a commensurate wavevec- 
tor Q = (±7r, ±7r), and an incommensurate one with 
Q = <5i(±7r, ±7r). The calculation of the real space 
Fourier transform finally yields 



x' G (r) = ^(Q)exp(-^), 



(8) 



where Ko is the modified Bessel function, and F(Q) = 
cos(Q x r x ) cos(Q y r y ) . 

Having determined the Ni induced Cu spin polarization 
(sf), it is straightforward to investigate the 17 NMR 
lineshape, determined by the coupling of the / = | nu- 
clear spins 17 I; to the Cu electron spins s, with spatially 
varying mean value (sf). The hyperfine Hamiltonian is 



fthf 



Heinle ^2 Ci ± Sl 



17 



(9) 



where 7 n ,7 e are the gyromagnetic ratios for the 17 nu- 
cleus and the electron, respectively. The hyperfine cou- 
pling constants Cu is dominated by a nearest neighbor 
hyperfine coupling C w 3.3T//is p^Sj - However, it was 
recently argued that a next-nearest neighbor hyperfine 
coupling C w 0.25C is relevant for the explanation of the 
spin-lattice relaxation rate jig] in I^-zSr^CuO^ We 
will therefore also consider its effects on the 17 NMR 
line. Using a mean field description of this hyperfine cou- 
pling by replacing by (sf ) of Eq. ||, we finally obtain 
for the shift of the resonance at a given 17 site r/ 



E C MX'(r t -r,). 



(10) 
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Here, the sum over i runs over the Cu spin sites, coupled 
to the 17 nuclear spin, whereas the j-summation goes 
over all Ni-sites. Furthermore, the constant prefactor A 
is given by ^ n ^eJhCcurieB /(gfiB) 2 ■ Note, v x as given 
in Eq. |To| is the shift of the 17 resonance with respect 
to the case without Ni impurities. 

To obtain the 17 NMR line shape, we create a his- 
togram l (y) = ^2rS(v — i>{) counting the number of 
nuclei with shift v. Since we want to compare the re- 
sulting distribution with the experimental data where 
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v-v (in arb. units) 

FIG. 2. The 17 lineshape l C onv{v)/Iconv{v = 0) for 
x — 2% Ni doping, C' — and two different values for the 
correlation length £: £ = 2 (dashed line), and £ = 4 (solid 
line). For both convolutions we used a = 20kHz 
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FIG. 3. Av(£)/Avo for a; = 2% Ni concentration and 
C = 0. The open diamonds present the results using ~x!q{(\), 
whereas the filled squares show the results obtained with 
Xl(q)- 



the line has a finite width even in the absence of im- 
purities, we convolute I {v) with a Gaussian distri- 
bution exp (— v 2 /(2a 2 )) /V27rcr 2 , yielding the lineshape 
I(v). By comparison with the experiments of Ref. ||] 
we expect that A^o = -y/21og 2a should be of the or- 
der of the high temperature (i.e. £ < 1) Ni-impurity 
induced linewidth. In the following calculation we there- 
fore choose a — 20kHz for both the Lorentzian and Gaus- 
sian x'(q). Finally, we define the resulting Av by half 
the width of the peak at half maximum. In Fig. we 
present the lineshape of the 17 O NMR signal, calculated 
with the Lorentzian form Xl(i) f°r two different values 
of £. We clearly observe that the line becomes broader 
as we increase £. From a comparison of Eq.(|lO|) with 
the experimentally measured broadening we can extract 
the value of the interaction J in Eq.(||). For C = 
and £(200K) = 4 (£(200K) = 3) we obtain J w 25 meV 
(43 meV). These values are accompanied by some uncer- 
tainties, but enable us to estimate the effects of a Cu spin 
mediated Ni-Ni spin (RKKY-type) interaction. We find 
within a self consistent mean field treatment of this in- 
teraction that the effect of the Ni-Ni interaction changes 
Av only within a few percent, consistent with the fact 
that no significant deviation from a Curie law was found 
in susceptibility measurements PJiTj]. 

In Fig. H| we present a comparison of Av(^) for the 
Gaussian x^(q) (open diamonds) and the Lorentzian 
Xl(i) (filled squares). Here we chose a Ni concentra- 
tion of x = 2%, C = 0, and Q = <5i(±7r,±7r) to be 
incommensurate with Si — 0.94 Q. We also compute Av 
for the commensurate case, and find that in general Av 
decreases. However, since the incommensurability, 1 — 8{, 
in YBa2Cu306.6, if present at all, is rather small, dif- 
ferences are negligible for £ < 8. In Fig. g we clearly 



see that the effect of Xc(q) and X^(q) on the behav- 
ior of the line width is qualitatively different. In agree- 
ment with the results by Bobroff et al. we find using 
Xc(q) that Av is basically independent of £ for all phys- 
ically reasonable values 2 < £ < 5. The Lorentzian form 
XL(q)j however, yields a much stronger increase in Av 
between £ = 2 and £ = 5 than the Gaussian. This result 
immediately implies that a temperature dependent £ is 
clearly compatible with the experimental results by Bo- 
broff et al. Furthermore, we find that the function /(£) 
of Eq. |l| behaves like /(£) ~ £ 3 / 2 in the Lorentzian case 
and /(£) ~ const, in the Gaussian case. This qualita- 
tively different behavior of Az/(£) for Xc( t l) anc l Xl( ( i) 
makes this experiment extremely sensitive to details of 
the momentum dependence of x'(q). 

Next we discuss the £ dependence of Av for different 
values of the Ni concentration x. We present our results 
for a Ni concentration of x = 0.5%, 2% and 4% and for 
C = 0.25C in Fig. |[ In agreement with the experimen- 
tal results we find that Av for a given £ increases with x. 
We believe that the results in Figs. || and [| also provide 
an explanation for the different behavior of Av in the un- 
derdoped (YBa 2 Cu306.6) and overdoped (YBa2Cu3C-7) 
samples. Bobroff et al. obtained that for the overdoped 
sample the variation of Av with T is much weaker than 
for the underdoped sample. As far as x'Cl) is concerned, 
the main difference between these two regimes consist in 
the value of £, namely £ = 1..2 for the overdoped and 
£ = 2. .4 for the underdoped sample. We see from Figs. || 
and [|that the £ variation of Av for the overdoped sample 
is much weaker than for the underdoped one, in agree- 
ment with the experimental results. 

Finally, we can use our numerical results to investigate 
in more detail the consequences of the T-dependence of 
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spin dynamics in the doped cuprates. A more correct 
description is provided by a Lorentzian-type form Xl(i)> 
which is fully compatible with the experimental data and 
a temperature dependent £. Though the resolution of the 
experiment does not allow us yet to determine the precise 
T-dependence of £, our analysis shows that £ consider- 
ably decreases with increasing temperature. 
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FIG. 4. Av(£)/Av for C = 0.25(7 and three different Ni 
concentrations: x — 0.5% (squares), x — 0.5% (triangles) and 
x = 0.5% (diamonds). 

TAv imp T 2G shown in Fig. g. Using g(£) cx £ |p) and 
the above results for /(£), it follows from Eq. yjfor the 
Gaussian case (/(£) ~ const.) that TAis unp T2G cx f , 
i.e. £ has to increase with increasing T. This result 
seems to be unphysical and thus strongly suggests that 
the Gaussian form Xcj( < l) ' s inappropriate for the de- 
scription of the spin susceptibility. On the other hand, 
for the Lorentzian case, /(£) ~ £ 3 / 2 and it follows 
TA^i mp T 2 G oc £ 1//2 , i.e. £ decreases as T increases, as 
we would expect. One can also solve Eqs.(|l|) and &j to 
obtain a as a function of T. However, our result possesses 
error bars which are quite large. The conclusion that a 
is independent of T is acceptable within those errorbars, 
however, a weak T dependence cannot be excluded. 

It is important to contrast our findings with the obser- 
vations of INS experiments. In YBa2Cu306+,5, INS ob- 
serves a T-independent broad peak around (n, 7[1 , above 
T c , resembling a Gaussian form of x"(q, ^) Jj],@]- How- 
ever, strong indications for incommensurate peaks with 
Lorentzian like shape in YBa2Cu30g.6 j§ suggest that 
the broad structure around (tt, 7t) is only a superposition 
of incommensurate peaks. Its width is therefore dom- 
inated by the merely T-independent incommensuration 
instead of This is consistent with the recent analy- 
sis by Pines [O that the overall magnitude of \" (<L w ) 
in YBa2Cu306+<5, as obtained from NMR experiments, 
necessitates a considerable improvement of the experi- 
mental resolution of INS experiments to resolve the in- 
commensurate peaks in the normal state. 

In conclusion we obtain from the analysis of the 17 
NMR data by Bobroff et al. and the T 2G data by Taki- 
gawa that the correlation length £ must have a substan- 
tial temperature dependence. A detailed analysis shows 
that the Gaussian form x'g( < i) °f the s P m susceptibil- 
ity can be excluded as an appropriate description of the 
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